Draft version February 6, 2010 

Preprint typeset using LAT^X style cmulatcapj v. 11/10/09 



O 

(N 



43 

6 

in 

c3 



> 

On 
O 

(N 
O 
O 



X 



ACHROMATIZING A LIQUID-CRYSTAL SPECTROPOLARIMETER: RETARD ANCE VS STOKES-BASED 

CALIBRATION OF HI VIS. 

D. M. Harrington 1 , J.R. Kuhn 1 , C. Sennhauser 2 , E.J. Messersmith 1 & R.J. Thornton 3 

Institute for Astronomy, University of Hawaii, Honolulu-HI-96822 and 
Institute for Astronomy, ETH Zurich, 8093 Zurich, Switzerland 
Draft version February 6, 2010 

ABSTRACT 

Astronomical spectropolarimeters can be subject to many sources of systematic error which limit 
the precision and accuracy of the instrument. We present a calibration method for observing high- 
resolution polarized spectra using chromatic liquid-crystal variable retarders (LCVRs). These LCVRs 
allow for polarimetric modulation of the incident light without any moving optics at frequencies 
>10Hz. We demonstrate a calibration method using pure Stokes input states that enables an achrom- 
atization of the system. This Stokes-based deprojection method reproduces input polarization even 
though highly chromatic instrument effects exist. This process is first demonstrated in a laboratory 
spectropolarimeter where we characterize the LCVRs and show example deprojections. The process 
is then implemented the a newly upgraded HiVIS spectropolarimeter on the 3.67m AEOS telescope. 
The HiVIS spectropolarimeter has also been expanded to include broad-band full-Stokes spectropo- 
larimetry using achromatic wave-plates in addition to the tunable full-Stokes polarimetric mode using 
LCVRs. These two new polarimetric modes in combination with a new polarimetric calibration unit 
provide a much more sensitive polarimetric package with greatly reduced systematic error. 
Subject headings: techniques: polarimetric 



1. INTRODUCTION 

Stellar high resolution spectropolarimetry is an under- 
utilized and often complex technique, in part because po- 
larimetric signatures are often less than 1% and in many 
cases less than 0.1%. Nevertheless, when it can be reli- 
ably applied it is a powerful remote sensing diagnostic. 
In this small-signal regime, where high signal-to-noise is 
essential, there are important instrumental limits to the 
derived accuracy and precision. 

There are several methodologies used for making pre- 
cision measurements. The most obvious but practically 
most difficult method is to build a specialized instru- 
ment (and telescope) specifically optimized for polarime- 
try. One minimizes the number of optical elements, keeps 
all optical folds to near-normal incidence, minimizes scat- 
tered light and has some form of modulation or beam- 
swapping to remove systematic effects. 

A spectropolarimeter has many functional components 
essential for sensitive linear or circular polarization mea- 
surements. For example, polarizing beamsplitters, fixed 
and variable waveplates, and specialized synchronous de- 
tectors are common sources of polarimetry-specific sys- 
tematic errors. These errors include detector variations 
(pixel-to-pixel efficiency and field illumination effects), 
changes in illumination during the exposures comprising 
a measurement (telescope guiding error, beam wander 
from moving optics, atmospheric transparency variations 
and seeing) and any polarimetric or chromatic effects in- 
duced by the entire optical system. Some of these are 
easily correctable and some are not. Typically, one or 
two retarders are inserted in the beam to modulate the 
incident polarization for detection. These can be sev- 
eral types of achromatic retarders, Fresnel rhombs, liq- 
uid crystal variable retarders (either ferro-electric or neu- 
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matic) or piezo-elastic modulators. In some cases, the 
instrument, retarder or analyzer rotates to accomplish 
the modulation. Various choices of retardance have been 
used for various applications to optimize sensitivity to 
different types of incident polarization or for different ex- 
pected signatures. Since typical polarimetric detections 
at these small levels are often limited by instrument sys- 
tematic effects, careful designs must mitigate as many 
error sources as possible. 

Notable examples are the ESPaDOnS fiber-fed spec- 
trograph on the 3.6m Canada France Hawaii Telescope 
(CFHT) and it's copy Narval on the 2m Telescope 
Bernard Lyot (TBL) telescope (cf. Semel et al. 1993, Do- 
nati et al. 1999, Manset & Donati 2003). These instru- 
ments work at a resolution R^65,000 covering ~3700- 
10500A and are in active use by many (cf. MAPF0 
or MiMES3 )■ The two telescopes are equatorial and 
the polarimetric modulation is done immediately after 
a collimating lens and an atmospheric dispersion cor- 
rector (ADC). The instrument is fiber- fed and has a 
time- variable continuum polarization at the several per- 
cent level (ESPaDOnS Instrument Website 0). A similar 
continuum polarization is seen in the William Wehlau 
Spectropolarimeter unit with a similar design to the ES- 
PaDOnS package (Eversberg et al. 1998). The trans- 
missive optics (lens and ADC) currently induce a time- 
variable cross-talk of a few percent, but for most spec- 
tropolarimeters this is quite benign. Lower spectral res- 
olution instruments, such as HPOL and ISIS, typically 
have more stable polarimetric properties that have yet to 
be seen at higher spectral resolution (Wolff et al. 1996, 
ISIS Spectropolarimetry Manual Tinbergen & Rutten 

1 http://lamwws.oamp.fr/magics/mapp/FrontPage 

2 www.physics.queens.ca/ ~ wade/mimes/ 

3 www.cfht.hawaii.edu/Instruments/Spectroscopy/Espadons/ 
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19970). 

Ultimately it is calibration of the telescope and instru- 
ment that achieves the highest possible system polari- 
metric accuracy. Even though it is difficult to keep in- 
strumental polarization below 1%, the effort in build- 
ing a low polarization system isn't lost, as the ultimate 
calibrated polarization performance improves multiplica- 
tively. Thus, sensitive techniques which determine the 
Mueller matrix of the telescope/polarimeter (cf. Beck et 
al. 2005a & b, Kuhn et al. 1994, Joos et al. 2006, Patat 
& Romaneillo 2006, Tinbergen 2007) will achieve a lower 
overal system noise when the raw instrumental errors are 
minimized. 

High-resolution astronomical spectropolarimeters have 
common design elements. Typically, rotating achromatic 
retarders are placed before a calcitc-based dual-beam an- 
alyzer - a Wollaston prism in collimated space or a Savart 
plate at an image plane. ESPaDOnS and Narval use two 
half-wave and one quarter-wave Fresnel rhombs before 
a Wollaston prism. Two other dual-beam instruments, 
PEPSI on the 8.4m LBT at R-310,000 and HARPS on 
the ESO 3.6m telescope at R~115,000 are in various 
stages of construction (Strassmeier et al. 2003 & 2008, 
Snik et al. 2008). 

The HARPS polarimetric package was required to fit 
inside an already-existing space. This required the ana- 
lyzer to be a Foster prism with a cylindrical lens on one 
beam and a CaF2 prism for the second beam. The de- 
sign includes separate quarter- wave and half- wave super- 
achromatic plates for circular and linear polarization that 
cannot be used simultaneously. 

The PEPSI design is similar to ESPaDOnS in that 
a lens collimates the telescope focus and retarders with 
a Wollaston perform the modulation. Another lens in 
combination with an atmospheric dispersion corrector 
form an image on fibers which feed the spectrograph. 
Instead of Fresnel rhombs, a super-achromatic quarter- 
wave plate is chosen. Linear polarization sensitivity is 
achieved by physically rotating the entire polarization 
package. 

The High-resolution Visible and Infrared Spectrograph 
(HiVIS) we discuss here was an instrument initially built 
for spectroscopy but modified for linear spectropolarime- 
try (Thornton et al. 2003, Harrington et al. 2006, Har- 
rington & Kuhn 2008). A Savart plate and a rotating 
achromatic half- wave plate were installed at the entrance 
slit. One of the main complications of this system is 
the eight oblique reflections between the sky and the po- 
larization analyzer. Since the 3.67m Advanced Electro- 
Optical System (AEOS) telescope is altitude-azimuth 
and HiVIS is in a coude path, the many oblique reflec- 
tions cause pointing-variable cross-talk that can com- 
pletely swap linear and circular polarization from tele- 
scope input to output at some wavelengths and point- 
ings. Nevertheless, the instrument has been successfully 
used for sensitive linear spectropolarimetric studies of 
stars as well as other studies (cf. Harrington & Kuhn 
2007, 2009a, 2009b). 

Rapid modulation of the incident beam polarization, 
synchronous with the detector, is the one technique for 
removing time-dependent systematic effects. If this mod- 
ulation exceeds a kilohertz, then even atmospheric seeing 



errors can be minimized. This technique has been devel- 
oped by the solar community. For instance, the vari- 
ous incarnations of the ZIMPOL I and II solar imaging 
polarimeter have used piezo-elastic modulators or ferro- 
electric liquid crystals in combination with charge shuf- 
fling on a masked CCD to remove seeing induced sys- 
tematic errors (Povel 2001, Stenflo et al. 1997, Gandor- 
fer et al. 2004, Stenflo 2007). The Advanced Stokes 
Polarimeter (ASP) and La Palma Stokes Polarimeter 
(LPSP) are other notable examples (cf. Elmore et al. 
1992, Lites 1996, Martinez Pillet et al. 1999). This tech- 
nique has been adapted for night-time spectropolarimet- 
ric use at the Dominion Astrophysical Observatory us- 
ing ferro-electric liquid crystals and a fast-shuffling un- 
masked CCD. Though this instrument can only record 
a small spectral region at lower spectral resolution, fast 
modulation removes several systematic errors (Monin et 
al. in Prep). 

We present here a spectropolarimeter using retardance- 
varying liquid crystals. These liquid crystals add ad- 
ditional complexity as the retardance is chromatic and 
varies as a function of temperature. However, these chro- 
matic liquid crystals can be quite useful when considering 
their performance characteristics. The retardance varies 
without physical motion of the LCVRs. This removes 
systematic errors caused by moving an optical element 
during or between exposures. Another advantage is that 
this type of LCVR can switch from zero to half-wave re- 
tardance in several milliseconds. These retarders can be 
used in any type of fast-switching application. Another 
trade-off lies in detector cost and complexity. In cross- 
dispersed echelle spectropolarimeters, the detector area 
must be used efficiently. If one can encode polarimet- 
ric information in a minimum number of spectral orders 
and minimize the required number of CCD reads, one 
can optimize a system for faster, more efficient operation. 
These trade-offs motivate the use of LCVR systems. 

We show here that a Stokes-based deprojection method 
allows reconstruction of the incident polarization even 
when using LCVRs. In order to demonstrate this de- 
projection and improve the utility of HiVIS, we have im- 
plemented two different polarimetric modes. The first 
is a standard full-Stokes polarimetric mode using two 
rotating achromatic wave plates. The second is a fast- 
switching polarimetric mode using the LCVRs. With 
this new full- Stokes capability we were able to thor- 
oughly calibrate and compare the polarization properties 
of HiVIS and quantify the sources of systematic error and 
polarimetric cross-talk caused by the various optical ele- 
ments. 

In this paper we will demonstrate the calibration of 
the new retarders, performance characteristics of the up- 
graded HiVIS, and the Stokes-based deprojection rou- 
tines we have developed. First, the performance and 
calibration of the retarders using a laboratory spectropo- 
larimeter will be presented in section 2. The calibration 
of the LCVRs as well as the Stokes-based deprojection 
methods will be presented in section 3. The upgraded 
HiVIS spectropolarimeter properties will be presented 
in section 4. Section 5 will present the HiVIS spectro- 
graph polarization properties and section 6 will present 
the AEOS telescope polarization properties. 
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2. RETARDANCE FITTING AND OPTICAL 
CHARACTERIZATION 

In the Stokes formalism, the polarization state of light 
is denoted as a 4- vector: 



I 

Q 
u 

V 



(1) 



In this formalism, / represents the total intensity, Q 
and U the linearly polarized intensity along position an- 
gles 0° and 45° in the plane perpendicular to the light 
beam, and V the right-handed circularly polarized inten- 
sity. Note that according to this definition, linear polar- 
ization along angles 90° and 135° will be denoted as — Q 
and —U, respectively. Furthermore, Q and U and V are 
said to be orthogonal states of polarization. The degree 
of polarization can be defined as a ratio of polarized light 
to the total intensity of the beam: 



V /Q2 + U 2 + V 2 



(2) 



The intensity of unpolarized light will be represented 
in this paper as: 



I u = (1 - P)I 



(3) 
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The normalized Stokes parameters are denoted with 
lower case and arc defined as: 



(4) 



For details on polarization of light, see an excellent 
text: Collett 1992. 

The upgrade to full Stokes polarimetry on HiVIS re- 
quired the addition of an achromatic quarter- wave plate 
to the standard polarimctric mode and the installation 
of a completely new mode using two LCVRs. New cal- 
ibration optics were also installed to deliver broad-band 
linear and circular polarization to the HiVIS spectro- 
graph. We acquired and tested two 'compensated' LRC- 
200 IR-1 Mcadowlark LCVRs. These retarders were de- 
signed to deliver zero to over half-wave retardance over 
the 6500A to 9500A range. We also acquired a Boulder 
Vision Optic quarter-wave plate (QWP) and half-wave 
plate (HWP) designed to be achromatic in the 4000A to 
7000A range. Calibration polarizers consisted of a pair of 
Codixx wire-grid polarizers (1" round), a pair of Edmund 
High- Contrast polarizers, a pair of Edmund wire-grid po- 
larizers (50mm square) and a pair of simple sheet linear 
polarizers (2" round). We will denote these as Codixx, 
Edmund, wire-grid and sheet polarizers respectively. 

We describe below the characterization of these new 
optics using our laboratory spectropolarimetcr. The 
spectropolarimeter is essentially a collimated halogen 
light source with a pair of polarizers to generate and 
analyze polarized light. Light is collected in a fiber- fed 
Ocean Optics USB 200 spectrograph. The polarimetric 
efficiency for all the polarizer pairs used in this paper 
were measured with the laboratory spectropolarimeter. 
The sheet, Edmund, Codixx and wire-grid polarizer pairs 



were first measured both crossed (minimum transmis- 
sion) and parallel (maximum transmission) to character- 
ize their degree of polarization. The wavelengths showing 
>95% polarization for each polarizer pair are shown in 
Table 1. 



TABLE 1 
Optic Wavelength Ranges 



Name 


Wavelength (A) 


Sheet 


4200-8200 


Edmunds 


5200-10000 


Codixx 


5200-10000 


Wire-grid 


3800-10000 


LCVR 


6500-9500 


Wave- plate 


4000-7000 



This Table shows the name of each optic used in calibrations. The 
top section describes the polarizers used and the measured range 
where the measured polarization was >95%. The bottom section 
shows the design range for the retarders used. 

The sheet polarizer has a high degree of polarization 
from roughly 4800A to 8200A. The sheet polarizer was 
>80% polarizing from the near-UV limit of the spectro- 
graph at 3800A, above 95% in the 4200A to 8200A region 
and falling sharply with wavelength after 8200A. Both 
the Edmund and Codixx polarizers show high degree of 
polarization longward of 5200A. The wire-grid polarizers 
repeatably showed between 99.0% and 99.8% polariza- 
tion for the full spectral range of 3800A to lOOOOA. All 
of the measurements were reproducible to better than 
1%. Our four well-calibrated polarizer pairs are used to 
test the optics in the laboratory and in the HiVIS spec- 
trograph. In order to illustrate the repeatability and ro- 
bustness of the measurements, two independent data sets 
were taken with different polarizers. The sheet polarizers 
were used for short wavelengths and the Edmund polar- 
izers were used for long wavelengths. 

2.1. Achromatic Wave-Plate Characterization 

To calculate retardance, we take the ratio of trans- 
mission spectra with the retarder in between polarizers. 
The maximum transmission spectrum Jy is recorded with 
parallel fast axes of both polarizers and wave plate. This 
yields the transmission through all optics of the system. 
Next, the fast axes of the retarder and second polarizer 
are rotated 45° and 90° respectively. This spectrum is 
labeled I±, to denote the crossed polarizers. We then 
define the wavelength dependent transmission spectrum 
as 



T(A) 



i±W 



(5) 



Thus the transmission spectrum is always between 0% 
and 100% and is independent of the spectrum of the light 
source. Any polarization-independent absorption caused 
by the retarder or polarizers is included in both Ij_ and I|| 
and does not influence the transmission spectrum. Since 
the fast axis alignment is independent of wavelength and 
the polarizers are used only in their highly polarizing 
wavelength ranges, the transmission can be directly re- 
lated to the LCVR retardance. 
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4000 6000 8000 10000 4000 6000 8000 10000 

Wavelength in Angstroms Wavelength in Angstroms 



Fig. 1. — This Figure shows the measured transmission spectra T(A) calculated as in Equation \E\ for the QWP and HWP (dot-dash lines). 
The corresponding relative intensity spectra Li and Ix are shown as solid lines. The HWP is shown in the top two panels and QWP is 
shown in the bottom two panels. The left hand y-axis is relative intensity. The right hand y-axis shows the transmission ratio in percent. 
The relative intensities for the HWP are recorded through crossed Edmund polarizers on the top left panel and crossed sheet polarizers 
on the top right panel. The QWP intensities are recorded through parallel Edmund polarizers on the bottom left panel and parallel sheet 
polarizers on the bottom right panel. There are four individual measurements repeated at 0° , 90° , 180° and 270° rotation for the polarizer 
pairs giving rise to the many curves in each panel. 



To describe how polarized light propagates through 
any optical system, a Mueller matrix is constructed 
which describes how every incident state is transferred 
to an outgoing state. The Mueller matrix is a 4x4 set 
of transfer coefficients which multiplies the input Stokes 
vector to create the output Stokes vector: 



ML 



(G) 



If the Mueller matrix for a system is known, then one 
inverts the matrix and dcprojccts a set of measurements 
to recover the inputs: 



I 



M L, 



(7) 



The inverse Mueller matrix for the system we call the 
deprojection matrix. One can represent the individual 
Mueller matrix terms as describing how one incident 
state transfers to another: 



M 



( II QI UI VI \ 

IQ QQ UQ VQ 

IU QU UU VU 

\ IV QV UV VV ) 



(8) 



In our lab spectropolarimeter 9 measures the angle in 
the plane perpendicular to the light beam. If we define 
9 = 90° to be the direction of Stokes Q, we can write 
the standard Mueller matrix for a wave plate oriented at 



9 = 45° with its own transmission function T ret and a 
variable retardance 6 as: 



M ret (45° 



/ 1 o\ 

cos0 — sin^ 

1 

V sin cj) cos 4> J 



(9) 



The Mueller matrix of an ideal linear polarizer looks 
as follows: 



Mp z — -^T po i 



( 1 Tl 0\ 
Tl 10 

\ 0/ 



(10) 



Where the negative sign corresponds to 6 = 0° while 
the positive sign is 9 = 90°. Since the light of the cal- 
ibration lamp is unpolarized, Ii am p = Ti amp [1,0,0, 0] T , 
the Stokes vector of the light incident on the retarder can 
be written as 



_ 1 lT 
— —TiampTpoii [1,— 1,0, 0J 



(11) 



Our analyzer (the second polarizer) has a transmission 
function T po i2 and is a polarizer oriented vertically. With 
this Mueller matrix and input Stokes vector, we solve the 
transfer equation to recover the phase variations: 
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M po/ M ret I 



retain 



(12) 



If we note that \li am pTpoiiT ret T po i 2 = I\\ we write the 
output Stokes vector as: 



( l + cos{4>) \ 
— 1 — cos((j}) 





(13) 



This represents a purely polarized output beam with 
I = —Q. Since we are measuring this intensity through 
crossed polarizers, the output Stokes vector is simply Ij_ 
as denoted above. The only variable is a wavelength de- 
pendent retardance (</>). The measured output intensity 
becomes: 



/ X (A) = /||(A)(1 



(14) 



In terms of the transmission function defined in Equa- 
tion [5] we get an equation for the retardance: 



- l [T{\)} 



(15) 



This equation applied to the measured transmission 
curves gives us the absolute retardance variation with 
wavelength (</>(A)). 
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Fig. 2. — The average deviation in retardance (in degrees phase) 
calculated for both quarter and half wave plates through polar- 
izer sets. The equation <f> = cos~ 1 [T(A)] is implemented for each 
observation. An arbitrary offset has been applied for clarity. 



The measurements of I|| , Ij_ and the corresponding re- 
tardances are shown in the four panels of Figure [T] The 
HWP measurements are shown in the top two panels. 
The QWP measurements are shown in the bottom two 
panels. Both wave plates were rotated through 360 de- 
grees for both polarizer sets giving four independent mea- 
surements. In order to additionally test the repeatability 
of the mounting and aligning procedure, the QWP was 
unmounted and remounted for an additional set of mea- 
surements. This remounted measurement was indistin- 
guishable from the others. 



An ideal HWP oriented with the fast axis at 45° to 
the input polarization rotates the input state to the ori- 
entation of the final polarizer, giving 100% transmission 
through crossed polarizers. An ideal QWP would only 
give 50% transmission through both parallel and crossed 
polarizers. Thus, the bottom two panels have the right- 
hand y axis going to 50%. Since the transmission spec- 
tra of the wave plates are only useful when the polarizers 
have a high degree of polarization, these transmission 
spectra are only plotted for long wavelengths in the left- 
hand boxes and are only plotted for short wavelengths in 
the right-hand boxes. 

Since the QWP gives 50% transmission through both 
parallel and crossed polarizers, we decided to leave the 
polarizers parallel for the QWP measurement as it in- 
volves one less rotation of an optical element and gives 
less chance for systematic error. Thus the top two pan- 
els of Figure [1] show measured relative intensities either 
near or 25 while the bottom two panels show relative 
intensities near 12 or 25. One can also see the degree of 
polarization of the polarizers decreasing outside their de- 
signed wavelength ranges since the measured intensities 
show only small differences between crossed and parallel 
configurations. 

The corresponding retardance variation, computed as 
in Equation [T51 can be seen in Figure [5J Both QWP 
and HWP phase deviations measured through both sheet 
and Edmund polarizer pairs are shown. The HWP shows 
less than 10 degrees retardance deviation for the 4000A 
to 8000A region. The QWP is slightly less chromatic at 
longer wavelengths but more chromatic short of 5000A. 

3. ACHROMATIZING AN LCVR POLARIMETER 

In order to use the LCVRs in a spectropolarimeter, 
one must know how the retardance varies with applied 
voltage and wavelength. The transmission function is 
highly chromatic but independent of polarization. Fig- 
ure [3] shows the measured LCVR transmission function. 
A Varian Cary 5000 spectrophotometer was used to mea- 
sure the transmission of the two LCVRs we label A and 
B. The Cary spectrophotometer was calibrated for dark 
and unobscured transmission. The LCVR transmission 
was measured at three separate rotation angles. Since we 
measured the Cary light source to be highly polarized, 
and the LCVRs are highly birefringent and chromatic 
with zero applied voltage, any rotational modulation of 
the LCVRs could come from the polarization properties 
of the incident light, detector or of the LCVRs them- 
selves. 

There were significant differences between the two 
LCVRs at short wavelengths but there was no significant 
dependence on LCVR orientation. Both LCVRs trans- 
mitted nearly perfectly longward of 6000A. LCVR A had 
a significantly lower overall transmission around 4000A. 
Both LCVRs showed small-amplitude ripples that were 
orientation dependent. LCVR A transmits over 50% 
more light at very short wavelengths, but LCVR B trans- 
mits significantly more from 4000A to 6000A. 

Using the lab spectropolarimeter we measured the 
transmission spectra, I|| , with the LCVR fast axis aligned 
with both the front and rear polarizers. The transmission 
spectra, I±, was then measured for various LCVR volt- 
ages between crossed polarizers with the fast axis of the 
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4000 6000 8000 10000 
Wavelength in Angstroms 

Fig. 3. — The throughput (transmission) of the two liquid crystal 
variable retarders measured by the Varian Cray spectrophotome- 
ter. Each LCVR was measured at three different rotation angles 
and they are nearly indistinguishable on this plot. The Cray inter- 
nal light source is highly polarized at most wavelengths but only a 
mild 'ripple' is seen at the different rotation angles. 



LCVR rotated by 45°. The retardance values were recov- 
ered using Equation[T5]for each voltage in the same man- 
ner as the achromatic wave plates. It should be stressed 
that when the LCVR fast axis is aligned with the axis 
of either polarizer, there is no measurable effect with ap- 
plied voltage. The LCVRs show no measurable rotation 
of the fast axis with voltage or wavelength. 

These compensated Meadowlark LCVRs are designed 
to have zero retardance around 5-7 volts with over half- 
wave retardance at low voltages. With these specifica- 
tions, we convert the measured transmission spectra in 
to retardance spectra. We find the half and full trans- 
mission points through crossed polarizers as a function 
of wavelength and voltage then assign the retardance as 
cos- x (T). 

To illustrate this process, another set of measurements 
of both LCVR A and B was taken with 56 separate 
voltages from OV to 10V. Sampling of 0.25-0. 50V per 
step was performed around regions of rapid change with 
wider sampling at low and high voltages. These mea- 
surements were in agreement with the ones presented 
previously. Several hundred spectra were again recorded 
and normalized for each voltage and averaged to create 
a single precise transmission spectrum. Figure U shows 
these transmission spectra interpolated to a regular volt- 
age and wavelength grid. The contours for half and full 
transmission are overlaid representing 1/4, 1/2 and 3/4 
wave retardances. 

The structure of Figure 2] arises from the liquid crystal 
material properties and design of the LCVR. The liquid 
crystal itself is a layer of uni-axial birefringent material. 
The molecular axes are naturally aligned and the mate- 
rial is fixed between two rubbed glass surfaces to orient 
the material. On the inside of the glass surfaces is a 
thin transparent layer of indium-tin oxide that allows a 
voltage to be applied across the material. The material 
has maximal birefringence under low voltage as all the 
material is aligned parallel to the rubbed glass surfaces. 



As the voltage across the liquid crystal increases, the liq- 
uid crystal material in the center of the layer is torqued 
in response to the electric field. This central material 
aligns itself with the field while the material attached to 
the glass surfaces remains fixed and supplies a restoring 
torque. As voltage is increased the birefringence is re- 
duced as the total degree of molecular alignment in the 
liquid crystal is reduced. The molecules in the center 
become aligned with the electric field while those at the 
surface remain aligned with the surface. Thus the elec- 
tric field scrambles the bulk alignment of the material. 
Since these liquid crystals cannot sustain a DC voltage 
for a long time without damage, an AC square wave at 
2kHz is applied which simply alternates the sign of the 
field while preserving the amplitude. 




1 2 3 4 5 6 
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Fig. 4. — The LCVR-A transmission spectra through both the 
Edmund and sheet polarizers interpolated onto a regular voltage 
grid. Contours representing maximum transmission (and 1/2 wave 
retardance) as well as half-maximum transmission (and 1/4 & 3/4 
retardance) are overlaid. Red corresponds to maximum transmis- 
sion (1) and black corresponds to zero transmission. 

One can see in Figure 0] that there is a region between 
and 1 volts where there is essentially no change with 
voltage. This corresponds to the region where the electric 
field across the liquid crystal is insufficient to cause any 
molecular rotation in the material. There is a strong 
change seen from IV to 4V as the liquid crystal molecules 
begin to move. Above 4 volts, the material is mostly 
saturated and there is essentially no birefringence and 
very low chromatism. 

The two LCVRs have very similar retardance at a fixed 
voltage. However, two-decimal place precision in voltage 
is required to accurately reproduce retardances and an 
independent calibration must be done for each LCVR. 
These retardance measurements are repeatable to bet- 
ter than a few percent. The results of Figure [4] were 
reproduced even after a complete rebuild of the lab spec- 
tropolarimeter. 

There is a significant temperature dependence for the 
LCVRs. Retardance measurements were taken at 20°C 
and 40°C. The computed retardance when heated was 
subtracted from the retardance measured when cool to 
form the temperature coefficients (in radians per de- 
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gree temperature change) shown in Figure O There is 
a decrease in retardance with an increase of tempera- 
ture ranging from 0.02 to 0.08 radians per degree. This 
range of voltages corresponds to the region of maximal 
sensitivity to LCVR input voltage, only spanning 1.2-3V. 
There is generally low retardance variation with temper- 
ature at higher voltages and longer wavelengths as the 
physical retardance values are actually lower. There is 
generally higher temperature sensitivity at shorter wave- 
lengths and lower voltages as these settings correspond 
to higher birefringence and stronger implicit voltage de- 
pendence. This shows that temperature is an important 
variable that must be controlled when using LCVRs. 
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Fig. 5. — The liquid crystal temperature dependence shown in 
radians induced phase change per degree temperature increase. All 
values are negative showing that retardance drops with increasing 
temperature. 



polarized beams which are sensitive to other Stokes pa- 
rameters in the incident beam. The typical setup with a 
Savart plate analyzer, shown in Figure [6j is to have the 
first LCVR with its fast axis aligned with the Savart plate 
axis and a retardance (j>\. The second LCVR is rotated 
by 45° and has a retardance 4>2- The typical procedure 
is to use combinations 0, 90, 180 and 270 degrees phase 
to achieve the modulation. 

The first step in this example is calibrating the liquid 
crystals. From measurements like those of Figure 0] for 
7590A one finds that voltages of 6.50, 2.80, 2.10 and 
1.65 correspond to the retardances of 0, 90, 180 and 270 
degrees phase respectively. These voltages then become 
the settings for a normal sequence. The LCVRs are set to 
their respective phases to measure one particular Stokes 
parameter, say (</>i, 02)=(O,9O) for +U and the exposure 
is taken recording the corresponding spectra. The LCVR 
phases for each setting can be denoted as: 
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This 



and 



. be more compactly written as 

02 where a and j3 denote the retardance angle at 
the nominal wavelength. The phases are denoted as 
a=(0,0,90,90,0,0) and /3=(0,180,90,270,90,270). Each or- 
thogonally polarized beam exiting the Savart plate is a 
pure polarization state represented as a scalar intensity 
coefficient multiplying a Stokes vector for pure +Q: 



3.1. Observing with an LCVR Spectropolarimeter 

As a way of illustrating the deprojection process we will 
illustrate the extraction of the Stokes parameters from a 
typical observing sequence for a hypothetical spectropo- 
larimetric observation of the 7590A TiO band. In a typi- 
cal dual beam spectropolarimeter sequence, one exposure 
will correspond to a measurement of one Stokes parame- 
ter. There are other modulation schemes that record lin- 
ear combinations of the Stokes parameters in other ways, 
but many night-time high-resolution spectropolarimeters 
function in this manner. 

In a typical sequence, a series of 6 exposures is taken 
corresponding to positive and negative Stokes parame- 
ters. The dual-beam analyzer gives two orthogonally po- 
larized spectra per exposure. The orthogonally polarized 
exit beams become the top and bottom spectra imaged 
on the detector. Since the analyzer sends / + Q to the 
top spectra and I— Q to the bottom spectra and all other 
states are split equally between top and bottom spectra, 
one can simply subtract the top spectra from the bottom 
spectrum to obtain a polarized spectrum in one Stokes 
parameter. 

The retarders act to modulate the incoming polariza- 
tion so that the analyzer will produce two orthogonally 
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What we are interested in is how the incident polariza- 
tion state maps to the intensity coefficients It op and lb ot . 
This can be calculated by simply multiplying the Stokes 
vector of the incident light by the combined Mueller ma- 
trix for the two LCVRs at their respective orientations 
and retardances and a linear polarizer which represents 
the Savart plate analyzer. The main result is that un- 
polarized light is always recorded with equal intensity in 
both top and bottom beams. For every LCVR retardance 
setting of the normal sequence, the desired Stokes param- 
eter is always recorded with a coefficient of ±1 while the 
other Stokes parameters are recorded equally with co- 
efficients of 0.5. The prescription for where each Stokes 
parameter falls in each spectrum and the two correspond- 
ing LCVR phases (shown parenthetically in degrees) is 
given by: 



Harrington, Kuhn, Sennhauser, Messersmith & Thornton 



LCVR1 LCVR2 



Savart Plate 



Fast 
Axis 



Fast 
Axis 



eo A 



oe 



Fig. 6. — The schematic for the Hi VIS spectropolarimeter. Light 
incident from the left hits the first LCVR, oriented with the fast 
axis vertical. The second LCVR is rotated by 45° . The Savart plate 
is two calcite crystals bonded such that the extraordinary beam for 
the first crystal becomes the ordinary beam in the second crystal. 
The Savart plate axis is oriented vertically, giving orthogonally 
polarized ±Q exit beams. 
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This 'normal sequence' records one Stokes parameter in 
one exposure by encoding that parameter in a single top 
or bottom spectra with no sensitivity to the other Stokes 
parameters. For instance, Stokes Q is entirely contained 
in the first exposure (htop and hbot) at the optimum 
wavelength and the other states are equally recorded in 
both beams. The htop spectrum minus the hbot spec- 
trum gives only Stokes Q sensitivity. The sign of each 
state is reversed in even versus odd exposure numbers. 

This sign reversal is the key to the dual beam concept 
as it enables systematic error subtraction and indepen- 
dence from detector effects. Each Stokes parameter is 
recorded on the same set of pixels with opposite sign. 
Subtraction of the two spectra removes detector noise but 
adds the Stokes parameters. Once this normal sequence 
is observed, we have recorded 12 spectra in 6 exposures 
where simple subtraction of individual spectra is enough 
to disentangle each Stokes parameter. The three expo- 
sures recording the negative Stokes parameters are used 
to remove CCD cosmetic effects and any pixel-to-pixel 
sensitivity variations. 

Even though the LCVRs are chromatic, the polariza- 
tion information at all other wavelengths is still present. 
Extracting the polarization spectra is possible but re- 
quires a more complex deprojection process. This gives 
rise to wavelength-dependent sensitivities for each Stokes 
parameter. In order to deproject the observations, mea- 
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Fig. 7. — The retardance versus wavelength for four selected 
LCVR voltages. The solid line shows LCVR A and the dashed 
line shows LCVR B. The vertical line shows the TiO wavelength 
of 7590A where the LCVR retardances were optimized to be 0.00, 
0.25, 0.50 and 0.75 waves. 



surements of the LCVR retardance must be used. Figure 
[7] shows the retardance for the four voltages used in this 
sequence. There is minimal difference between individ- 
ual LCVRs as the solid and dashed lines overlay quite 
nicely. In order to represent our two LCVR polarimeter 
for all wavelengths, we calculate the Mueller matrix for 
a 2-LCVR system with arbitrary phase. We will denote 
the cosine of a phase as c0 and the sine of a phase as s0. 
The Mueller matrix of the two LCVR system (M-2LG'vr) 
is just M ret (45,0i) multiplied by M ret (O,02): 
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(18) 

This can be multiplied out to create the full Mueller 
matrix for the liquid crystals: 
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The dual-beam spectropolarimeter essentially has two 
separate analyzers. The top and bottom spectra can be 
represented as having passed through a linear polarizer 
oriented vertical (0°) and a second linear polarizer ori- 
ented horizontal (90°). The spectra recorded through 
each beam can be calculated by multiplying the 2-LCVR 
Mueller matrix by the Mueller matrix for the linear polar- 
izers which select the individual top and bottom beams. 
If the beam recorded through the vertical polarizer is the 
'top' spectra and the horizontal polarization is the 'bot- 
tom' spectra, the equation for the intensity coefficients 
of the recorded spectra becomes: 



hop \ _ 1 
hot 



4 - C02<5 - S02S0lC/ + S0 2 C0lV 

4 + c(f} 2 Q + s(j) 2 s(t>iU - s0 2 c0iV 



( 2 °) 

With the equations for the coefficients multiplying each 
Stokes parameter in each of the exposures, typically 
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called the efficiencies, one can determine the chromatic 
effects of the LCVRs and how to deproject each of the 
Stokes parameters. Figure [8] illustrates the wavelength 
dependence of these efficiencies for each Stokes parame- 
ter optimized for 7590A. In the normal sequence, the top 
and bottom spectra record one Stokes parameter multi- 
plied by an efficiency ±1. The other Stokes parameters 
are all equally present in both top and bottom spec- 
tra with a coefficient of 0.5. Each panel of the Figure 
shows the efficiencies for an individual Stokes parameter 
as a function of wavelength. During the first exposure, 
shown in the top box, Stokes Q is correctly multiplied 
±1 at 7590 A as both LCVRs are at zero retardance and 
show low chromatism at high voltage. The measurement 
of — Q (dashed curve) puts the second LCVR at half- 
wave retardance and lower voltage. The chromatic ef- 
fects are strong enough to make both efficiencies +1 by 
4500A. Stronger chromatic effects are seen in the other 
Stokes parameters because of the general higher retar- 
dances (and lower voltages) present on both LCVRs. 




5000 6000 7000 8000 9000 
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Fig. 8. — The LCVR efficiencies from Equation 1201 as functions 
of wavelength using the first three retardance settings of the 'nor- 
mal' sequence. The top box shows the standard 'Q' sequence, the 
middle box the '[/' sequence and the bottom box the 'V sequence. 
The solid lines are the + Stokes parameters while the dashed lines 
are the - Stokes parameters. Only for a narrow wavelength range 
around 7590A does the normal sequence properly map incident 
polarization to the required ±QUV as noted by the vertical lines. 



3.2. Retardance Fitting Deprojection 

Given that we have now measured the retardance for 
the LCVRs for all possible applied voltages, an easy de- 
projection to attempt is simply to calculate the Mueller 
matrix of the LCVR polarimcter for all wavelengths us- 
ing Equation [TH] We will investigate a sample observing 
sequence for input linear polarization as an example of 
this deprojection via retardance fitting. 

In order to properly deproject a series of 12 spectra 
from 6 exposures into a single polarization measurement, 
one must create a 'deprojection matrix'. Typically, one 
calculates the retardance of the LCVRs and then solves 
the system of equations which represent the polarized 
radiative transfer for the spectropolarimeter. We will il- 
lustrate this process for 7590A-optimized measurements. 
The observing sequence gives a set of redundant equa- 
tions for the polarized spectra: 



These equations depend on the scalar intensity co- 
efficients and can be expressed in terms of sines and 
cosines of the LCVR retardances. The retardances for 
each LCVR (<j>i,(j>2) were expressed in Equation [TBI To 
simplify the notation, we will number the spectra by ex- 
posure (1 to 6) and note their detector location (top or 
bot). We note that the retardances are implicitly func- 
tions of wavelength. The coefficients are: 
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( 22 ) 

In order to solve this system we simply invert the equa- 
tions: 
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This is easily implemented for all wavelengths using 
the retardances we have calculated in the previous sec- 
tion. An example of the deprojection process using the 
lab spectropolarimeter will be shown here. The LCVR 
fast axes were positioned as a polarimeter with LCVR 
A at 0° and LCVR B at 45°. Eight complete linear 
polarization states were input by rotating the front po- 
larizer through 360° in 45° increments. A full polari- 
metric sequence was taken for each input state using the 
voltages chosen to optimize the polarimeter for 7590A. 
The overall transmission used to normalize each measure- 
ment was recorded through parallel polarizers with both 
LCVRs aligned with the polarizers. Since the retardance 
with wavelength functions are known, the complete po- 
larimetric properties of the input can be deprojected us- 
ing Equation [23] Figure [9] shows the deprojection for all 
eight input states. The input linear polarization was ex- 
tracted well for all wavelengths. The input states were 
recovered and the cross-talk states were generally around 
0.1 or less though with some chromatism and some values 
as high as 0.2. 

One can see from this simple example of a normal ob- 
servation sequence that there are still significant system- 
atic errors present from phase fitting. The LCVR phases 
are measured in a different optical setup than when the 
LCVRs are used as a polarimeter. This leads to a wide 
range of potential errors. As we shall see in the next sec- 
tion, a method where phase fitting is not required per- 
forms much better. 

3.3. Stokes-Based Deprojection 

A deprojection method which uses measurements of 
known input states is straightforward and direct. It in- 
corporates all effects of the optical path and not just 
those of the LCVRs. There is no need to calculate the 
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phase of the LCVRs and no interpolation is necessary. 
Using our lab spectropolarimeter, this is illustrated in 
the following example. 
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Fig. 9. — The deprojection of pure lin ear polarization input using 
the standard deprojection of Equation 1231 Though chromatic, the 
deprojected outputs shown here reproduce the pure input states 
with values near 1 . The cross-talk Stokes parameters are the curves 
that are below 0.3. The discontinuity around 7500A is the wave- 
length transition between sheet and Edmund polarizers. 



optical systems this definition is the most useful for sys- 
tem optimization and error propagation calculations (cf. 
Complain et al. 1999, DeMartino et al. 2003). The 
condition number is well below 10, and increasing only 
for wavelengths short-ward of 6000A, where the polar- 
ization performance in the system starts to drop. As two 
examples, Figure HU1 shows the Stokes-based deprojection 
method with the front polarizer input states orthogonal 
to the calibrations (— Q and —U respectively). The re- 
covered Stokes parameters are the solid top lines near 
1.0. The dotted and dashed lines show linear and circu- 
lar cross-talk. 

The circular cross-talk is near zero with the linear 
cross-talk value being around 0.05. This represents 
roughly a factor of 4 improvement from the method uti- 
lizing LCVR phases presented previously. It should be 
pointed out that the linear cross-talk is essentially achro- 
matic. This can be easily explained as simply an error in 
the rotation angle of the linear polarizer when set a neg- 
ative input state. With the input state being reproduced 
near 1.0 and the cross-talk states being 0.05 or less, the 
angular error computed as 0.5 tan -1 (input / crosstalk) 
is less than 1.5°. This corresponds to the error what one 
would expect for a manual rotation of the linear polar- 
izer. The simple fact that the cross-talk is achromatic 
suggests that it is the improper rotation angle of the po- 
larizer used to create the pure input Stokes parameters 
that is the cause of the error, not the improper removal 
of chromatic LCVR retardance. 



We inserted fully polarized light into our lab spec- 
tropolarimeter as pure inputs +Q, +U, +V (as well as 
— Q, — U, — V to show consistency). For each input state, 
we took the six measurements applying the 'normal se- 
quence' of voltages to our two LCVRs. In the lab 
spectropolarimeter we have a polarizer as an analyzer 
so we only record the 'top' spectra from Equation [22l 
Since our input space is three-dimensional, our output 
(prime) space will have (at most) three independent vec- 
tors, which we created by defining Q' = l2top — htop, 
U' = htop - htop, and V = htop - htop- Note that we 
could have chosen completely different combinations of 
measurements to build up our three prime space vectors, 
but our choice was influenced by our knowing that at 
7590A, these three vectors represent the (unnormalized) 
values of Q, U , and V from our retardance measure- 
ments. As orthonormal input basis vectors we use fully 
polarized light - all three pure Stokes parameters define 
the basis vectors. The two front linear polarizer states 
defining +Q and +U were used as well as a quarter wave 
plate oriented to feed in +V. These are our pure input 
states used for the Stokes-based deprojection. Then the 
columns of the Mueller matrix for our setup are directly 
given by (Q',U',V ) for each of the input basis vectors. If 
we now want to retrieve the input state of any measure- 
ment we again build up the quantities htop- htop, htop- 
htop, htop- htop, and we apply the inverse of our Mueller 
matrix to them. The problem turns out to be always 
well-posed, as one can see by looking at the (wavelength 
dependent) condition number of the Mueller matrix. The 
condition number, defined as the absolute value ratio of 
the highest to lowest eigenvalues of the mueller matrix, 
is a measure of the invertability of a linear system. In 
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Fig. 10. — The deprojected Stokes parameters are shown here. 
The input state is reproduced at a value near 1. The cross-linear 
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3.4. Benefits of Stokes-Based Deprojection 

The main benefit of Stokes-based deprojection can be 
seen when considering how errors effect measured Stokes 
parameters. In dual-beam spectropolarimeters, the stan- 
dard observing sequence is followed in order to accom- 
plish two independent goals. One goal is to create po- 
larization spectra that are independent of pixel-to-pixel 
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efficiencies. Since these pixel-to-pixel variations are usu- 
ally not known to better than a percent, precision spec- 
tropolarimetry requires observations to be independent 
of the gain. The other main goal is to choose retardances 
that map each Stokes parameter into these normalized 
double differences. For the standard sequence used with 
the achromatic wave plates this can be implemented as 
follows: 

Q 1 / Iltop — Iltop Ilbot — l2bot s 1/ x 

q = T = 9(7 777 7 777 — ) = 9 (?1 + q2 > 

One can show that this method is independent of pixcl- 
to-pixel gains by propagating gain coefficients through 
Equation 1241 Suppose that the incident radiation inten- 
sity in a top pixel is multiplied by a gain factor g t . On a 
bottom pixel the gain is gf,. In the standard choice of re- 
tardances, the intensity of the incident radiation switches 
beams so that the l\top intensity before detection is the 
same as l2bot- L&oi is identical to l2t op . Thus there are 
only two incident intensities denoted as It and lb which 
switch places between exposure. The incident radiation 
is modulatd so that the gains cancel and the normalized 
double difference (Equation [M]) becomes: 

q _ 1 , 9th ~ 9th 9bh - 9bh , _ h - h ^ 
2 g t I t + g t h 9bh + 9bh h + h 

This is exactly as desired - the pixel-to-pixel gain vari- 
ations have been eliminated and all that remains is a nor- 
malized difference which records one Stokes parameter. 
There is no systematic-error induced by pixel-dependent 
variations. The difference ratio method has the added 
benefit of removing the optical mis-alignment errors. The 
slit tilt and 'derivative error' caused by wavelength mis- 
alignments between 'top' and 'bottom' spectra are re- 
moved because the same pixels are used for each normal- 
ized difference (qi and q2). 

While the normalized double differences always can- 
cel the pixel-to-pixel gains, the chromatic effects of the 
LCVRs require a deprojection from input to measured 
QUV. The Stokes-based deprojection utilizes calibra- 
tion observations taken with as similar an optical con- 
figuration to the measurements as possible. One can see 
the quality of the Stokes-based deprojection by compar- 
ing the retardance-fitting deprojection of Figure [9] with 
Stokes-based deprojection of Figure [TO] Aside from the 
small ripples inherent in the LCVR transmission func- 
tion, the Stokes-based method is much more achromatic, 
less prone to optical sources of error and easier to imple- 
ment. 

4. HIVIS LCVR SPECTROPOLARIMETER 

After demonstrating Stokes-based deprojection with 
the lab spectropolarimeter, we installed an LCVR spec- 
tropolarimeter unit and the associated calibration optics 
on the HiVIS spectrograph. We are presently upgrad- 
ing HiVIS to include better optics, optomechanics and 
a fast-switching mode utilizing these LCVRs. This up- 
grade involves using a new charge- shifting detector. The 
original detectors, two 2k by 4k CCID-20 were replaced 
with new CCID-20s and electronics that allow charge to 
shuffle in one direction. If the modulation is performed 



faster than telescope guiding errors, a major source of 
systematic error can be removed. 

In order to do this fast modulation, the charge will be 
shifted on the CCD without reading out an exposure. 
This charge shifting must also be synchronized with the 
polarimetric modulation. The spectral orders will not be 
allowed to cross in the charge shifting process as there 
is presently no shuttering during an exposure. In order 
to maximize signal-to-noise and make efficient use of the 
telescope, the number of imaged spectra must be mini- 
mized and this optimization procedure implemented. We 
do not need to have one exposure per Stokes parameter. 
As long as the polarization information is encoded in the 
recorded spectra in an extractable manner, the number 
of exposures can be reduced. 

The new HiVIS spectropolarimeter has several up- 
grades. The additions include a second rotating Boulder- 
Vision-Optic (BVO) achromatic quarter- wave plate, are- 
mounting of all the polarizing optics with full xyz control, 
a new dekkar with full xyz control and polarization cali- 
bration optics. A new Savart plate was acquired with a 
larger clear aperture, more displacement and is more ac- 
curately bonded and a new mount which allows for easy 
switching between the standard wave plate mode and the 
fast-switching liquid-crystal mode. 

The new Savart plate has a 17% larger displacement 
and no detectable leak into the undesired e-e and 0-0 
beams (compared to the 1% measured for the old Savart 
plate shown in Harrington & Kuhn 2008). The new 
dekkar and Savart plate allows for more sharply defined 
and wider spectral orders on the focal plane. The dekkar 
stage allows for calibration optics to be mounted just up- 
stream of the slit and all pure input states (QUV) are 
reproduced with the achromatic wave-plates above the 
95% level for all wavelengths observed (5500A-7000A). 
Once this accurate reproduction of incident polarization 
at the slit was verified using the wave-plates the new 
2-LCVR spectropolarimeter fast switching system was 
installed. 

We first show the properties of the HiVIS liquid-crystal 
system and then illustrate the Stokes-based deprojection 
for the system. The variations in signal-to-noise ratio 
with wavelength in the deprojected output will be dis- 
cussed along with implications for the design and opti- 
mization of an achromatized spectropolarimeter. 

4.1. HiVIS LCVR Observations 

In order to make the HiVIS LCVR system very robust, 
empirical calibration hardware and software has been de- 
veloped. A new calibration unit consisting of an Edmund 
high-contrast polarizer and a BVO quarter-wave plate 
creates pure input states (±Q, ±J7, ±V) at the spectro- 
graph slit mirror. These pure inputs are used to derive 
the retardances and Mueller matrix coefficients of the 
LCVRs exactly as mounted and run during the night. 
As an example of this technique, we will use a sample 
set of H a observations and calibrations. These observa- 
tions were obtained using the Apogee 3k 2 detector. The 
Apogee has 13 spectral orders in the 5500A to 7000A 
range each with 3056 independent polarization measure- 
ments per order. We will define Stokes +Q as parallel to 
the slit. 

The LCVR voltages were selected to be only roughly 
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close to the 'normal sequence' for H Q . Deliberate round- 
ing errors in voltage were introduced to show that, even 
with significant retardance selection error, the Stokes- 
based deprojection method works well. Voltages of 7.25, 
2.95, 2.10 and 1.90 volts were chosen corresponding to 
retardances of -0.004, 0.215, 0.425 and 0.596 waves when 
mounted in the lab spectropolarimeter. This lab spec- 
tropolarimeter differs from the HiVIS setup in temper- 
ature and the two LCVRs were remounted introducing 
potential alignment variations. These values are rough 
and the deprojection does not require a derivation of the 
actual retardances. 
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Fig. 11. — This Figure shows the 6 measured polarization spec- 
tra of the normal sequence. The six individual ex pos ures produce 
12 spectra which are differenced as in Equation 1241 representing 
iQUV. The Stokes parameters have been normalized. The solid 
lines show ±q, the dashed line shows ±u and the dot-dash line 
shows ±1). To illustrate the repeatability of the measurement, two 
complete observation sequences were taken with pure +Q and — Q 
inputs. The sign has been reversed for all measured output states 
for the — Q input. Thus there are 24 spectra plotted here. See text 
for details. 



For each pure input state, the standard set of ±QUV 
observations listed in Equation 1161 were taken. The data 
were processed using the IDL reduction scripts described 
in Harrington & Kuhn 2008. We found that each echelle 
order required using a dark subtraction from the non- 
illuminated regions adjacent to each spectral order. A 
spatially and temporally variable dark current at the 
100- ADU level is present in the Apogee detector and is 
removed independently for each exposure. 

If the usual dual-beam polarization analysis is per- 
formed using these observations then the input pure po- 
larization states are poorly reproduced, as expected. The 
LCVR chromatism and our deliberate introduction of im- 
perfect retardances causes significant changes in how the 
Stokes parameters are recorded on the detector. This 
is the perfect situation to illustrate the power of this de- 
projection technique to correct an imperfect polarimeter. 
Figure fTTI shows the measured Stokes parameters derived 
from the six exposures for two full observing sequences. 
Pure +Q was input for one sequence and pure — Q was 
input for another sequence. The sign of the measured 
Stokes parameters has been reversed for the —Q input se- 
quence for ease of plotting. The high voltage setting cor- 
responding to zero retardance (7.25, 7.25) shows a near- 



perfect reproduction of ±q for all wavelengths (solid line 
on top) whereas the half-wave retardance on the second 
LCVR induces a fairly chromatic dependence with only 
70% reproduction of ±17. The u states are nearly zero 
at H a but are highly chromatic in that the dashed lines 
cross zero but reach -.4 and +.8 at short wavelengths. 
The v states are very similar to u but with a significant 
offset - there is no wavelength where both + and - v are 
zero. 

4.2. HiVIS LCVR Stokes-Based Deprojection 

Pure input states [+Q, +U, +V] were observed with a 
full polarization sequence described above (6 exposures, 
12 spectra). As each set of pure input states are ob- 
served twice, one set with lower retardances and one set 
with higher retardances, we actually get two independent 
measurements of each input state and hence two indepen- 
dent deprojection matrices. The measurements of each 
pure input state at low and high retardance form a trans- 
fer matrix at low and high retardance respectively. We 
then input pure negative states for testing the calibration 
routines [— Q, — U, — V]. These calibration observations 
also give us 6 exposures and 12 spectra which correspond 
to QUV measurements at both low and high LCVR re- 
tardance values. The transfer matrix calculated using the 
positive states were inverted and used to deproject each 
input negative state. The resulting empirical calibration, 
shown in Figure [T2l shows the residuals are less than 0.1. 
The three input states (negative pure states) map to es- 
sentially pure measured output states (the curves near 
+ 1). The sign was reversed for clarity. The curves near 
zero correspond to the circular 'cross-talk' that has not 
been removed by this deprojection procedure. 
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Fig. 12. — This Figure illustrates the empirical inversion process 
with LCVRs on HiVIS. Input states +Q, +U and +V were ob- 
served to create the deprojection matrix. The calibration optics 
were rotated to input negative states. These measured and depro- 
jected negative Stokes parameters are shown here. The input state 
is reproduced at a value near 1 with the cross-linear and circular 
states are quite small. 



These cross-talk values are quite small, being typically 
around 0.05 and show little chromatic variation. Consid- 
ering that the calibration optics were manually rotated 
and are on simple post-mounts, we consider this correc- 
tion to be excellent given the chromaticity of the LCVRs. 
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As noted earlier, achromatic offsets in cross-talk can re- 
sult from imprecise rotation of the calibration optics. It 
is the chromatic variation of the deprojectcd Stokes pa- 
rameters that is directly caused by errors in the depro- 
jection process. The chromatic variation in Figure [T2l is 
much smaller than the achromatic offsets from zero. 

4.3. HiVIS LCVR Optimization 

Since deprojection removes the constraint that one 
Stokes parameter correspond to one exposure, we are 
free to choose retardances that give the desired preci- 
sion for each Stokes parameter. One can optimize the 
observing sequence and retardances for a given require- 
ment by controlling the properties of the deprojection 
matrix. For instance, the condition number of the de- 
projection matrix gives the ratio of the maximum and 
minimum eigenvalues (cf. Compain ct al. 1999, De Mar- 
tino et al. 2003). This number is an upper limit to 
the maximum relative variance of the corresponding de- 
projected Stokes parameters. Said another way: Given 
equal exposure times for the measured spectra, the sta- 
tistical noise in the deprojected spectra will have a ratio 
of maximum to minimum variances that follows the con- 
dition number. We calculate the condition number by 
first applying the IDL single- value-decomposition routine 
SVDC to our deprojection matrix then taking the ratio 
of maximum to minimum eigenvalues (£ m m, £, m ax) ■ The 
condition number, c, is defined as: 



(26) 



We have investigated this relationship using our HiVIS 
observations. Using the deprojection matrix we have cal- 
culated for the LCVRs, we computed the condition num- 
bers as in Equation 1261 The condition number for the 
low-retardance setting (high voltages) was always below 
3. The condition number for the high-retardance set- 
ting (low voltages) was generally higher but always less 
than 4. This implies that, for a constant exposure time 
and statistical noise per recorded spectrum, the signal- 
to-noise of the resulting deprojected spectra will vary by 
no more than y/3 and y/Z respectively for the two set- 
tings. 

We then compared the variance of our polarimctric 
measurements both pre- and post- deprojection. Since 
the Apogee detector has 13 spectral orders each with 
3056 independent polarization measurements, we get 13 
variance measurements in the 5500A to 7000A range. 
The variance o~ pre ^ were computed for each order i before 
deprojection. Since the exposure times were constant, 
o~ P re,i = o~ pre . Accordingly, the variances for each order 
j post- deprojection <J p0 st,j were computed. Since the 
condition number gives the ratio of maximum to mini- 
mum variance, we define the maximum measured vari- 
ance change as 



Sv r , 



MAXji pP2fhl) 

(J p re 



(27) 



In addition, we also simulated error propagation 
through our calculations with simulated input Stokes pa- 
rameters. In the IDL reduction package, we simply re- 
placed the measured Stokes parameters (13 orders, 3056 



wavelengths per order) with noise about zero with a vari- 
ance of 1 (o~ pre =l). The deprojection was performed and 
the variances CpoTtj- were determined. We can compute 
the variances of this simulated data in analog to Equa- 
tion EH 



6v sim = MAX^Zj) (28) 

Figure [13] shows the condition number, measured vari- 
ance change (Sv meas ) and simulated variance change 
{5v S i m ) which track each other to better than 0.2 across 
the recorded wavelength range. 

The final test performed was to verify that the coeffi- 
cients in the deprojection matrices agree with the relative 
error values. The deprojection process multiplies the in- 
verse Mueller matrix (M -1 ) by the measured polarization 
spectra. Each deprojected spectrum is a sum of three in- 
put spectra multiplied by the corresponding deprojection 
matrix element. The variance of statistical noise will in- 
crease as the square of the deprojection matrix element 
and the variance will also increase as the sum of these 
squared values. The variance should go as n=a 2 +b 2 +c 2 
where a, b and c represent the row of the deprojection 
matrix (M _1 ). This we term the 'noise amplification co- 
efficient'. This coefficient also agrees with the simulated 
variance and the condition number to better than 0.2 as 
seen in Figure 1131 If we define each term of the depro- 
jection matrix as M" 1 for the i th row and ] th column, we 
can define these noise amplification coefficients n^ as: 



1\2 



(29) 



Since these noise amplification coefficients are a direct 
measure of the signal-to-noise computed for each depro- 
jected Stokes parameter, they can be used to optimize 
an observing sequence for any desired outcome. 
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Fig. 13. — This Figure illustrates the relationship between the 
condition number, measured variance, simulated variance, and 
'noise-amplification-coefficient'. Since the variances are computed 
as ratios, all four curves should agree. 



There are many schemes for tuning polarimetric mea- 
surements where the properties of the retarders (fast-axis 
orientation and retardance) are varied in order to give 
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some specified performance of the condition numbers. 
There are efficiency-balanced schemes, which maintain 
condition- numbers of 1, optimize for equal sensitivity be- 
tween linear (QU) and circular (V) or optimize on other 
constraints, (c.f. Nagaraju et al. 2007, Complain et al. 
1999, DeMartino et al. 2003). 

The ultimate limitations of this system lie in the mea- 
surement of the Mueller matrix terms. In this example, 
even where manual rotation of two different optics in- 
duced a substantial systematic error, the achromatized 
performance was acceptable. The measured cross-talk 
error in the output Stokes parameters was satisfacto- 
rily achromatic. Uncertainties in voltages or calculated 
LCVR phases do not influence the performance of the 
Stokes-based deprojection. There were systematic off- 
sets that can be explained as 1-2° errors in rotation of the 
calibration optics. The chromatic error in the recovered 
polarization are significantly lower. A motorized stage 
for the calibration optics should improve the achromati- 
zation even further. 

This system is designed to be used for line spectropo- 
larimetry which is an inherently differential technique. 
The continuum polarization accuracy of this system is 
limited by these deprojection uncertainties, as well as 
the polarization induced by the optics of the system. 
However, this calibration is quite sufficient to create a 
stable polarimetric reference frame for measurements of 
polarization across individual spectral lines. One of the 
most widely used high-resolution spectropolarimeters, 
ESPaDOnS, has a cross-talk of 2% and a wavelength- 
dependent, time- variable continuum polarization (Do- 
natietal. 2006, ESPaDOnS Instrument Website 0). Ab- 
solute polarization measurements suffer from many sys- 
tematic errors but the change in instrument properties 
over an individual spectral line are negligible. Hence, 
achromatization of the instrument performance to an ac- 
curacy of a few percent is quite sufficient to allow for 
accurate differential measurements. 

5. POLARIZATION PROPERTIES OF THE HIVIS 
SPECTROGRAPH 

In order to perform accurate polarimetry of astronomi- 
cal sources, a polarization calibration of the telescope and 
spectrograph optics is also required. The deprojection 
procedure we have just outlined is quite general and can 
also be applied to undo cross-talk induced by other opti- 
cal elements in an instrument. The AEOS telescope with 
5 oblique fold mirrors and the HiVIS spectrograph with 
its image rotator is a suitable example. The AEOS tele- 
scope induces pointing-dependent cross talk. The image 
rotator in HiVIS induces rotation-dependent cross-talk. 
The spectrograph itself also induces cross-talk even with- 
out the image rotator. Since the Savart plate analyzer is 
just after the slit, only the optics before the slit induce 
polarimetric effects. In this section, we will use the achro- 
matic wave-plate spectropolarimeter to measure the po- 
larization properties of the HiVIS spectrograph optics. 

The image rotator mount was aligned and remounted 
to allow repeatable removal and re-insertion of the im- 
age rotator with negligible influence on the spectrograph 
beam path. The optical path from the last coude mirror 
through the spectrograph fore-optics to the analyzer has 

5 www.cfht . hawaii.edu /Instruments / Spectroscopy /Espadons / 



many elements. There are three near-normal incidence 
fold mirrors and a 5.5m focal length collimator before 
coming to the image rotator (3 oblique reflections). An- 
other fold, the tip-tilt mirror, a 1.1m focal length sphere 
and another fold are used to create the stellar image at 
the slit. To calibrate with pure linear polarization, the 
wire grid polarizer was mounted at the calibration stage 
inside a mask to polarize the diffusion screen and to block 
unpolarized light. This serves as a polarized flat field cre- 
ating pure QU inputs to the spectrograph fore-optics. 



1.0 
0.9 
0.8 
0.7 
0.6 
0.5 

0.2 

0.1 

0.0 
-0.1 
-0.2 

1.0 

0.5 

0.0 
-0.5 
-1.0 
5500 



U Input 



Q Input 



6000 6500 
Wavelength in Angsroms 



7000 



Fig. 14. — This Figure shows the output Stokes parameters for 
four inputs (+Q, +U, —Q, —U) with the image rotator in the nom- 
inal orientation. Solid lines are for ±Q inputs, dashed lines are for 
±C/ inputs. The top box shows the measured output state for the 
same input state (eg. +Q out for -\-Q in). The middle panel shows 
the non-input linear Stokes parameter showing linear cross-talk of 
up to 0.2. The bottom panel shows circular polarization and is 
entirely cross-talk. 



5.1. Polarization With Image Rotator 

The polarization response of the spectrograph fore- 
optics was first measured with the image rotator in the 
beam and aligned vertically. Pure ±Q and ±U states 
were input. The cross-talk effects can be best illustrated 
by separating the output measurements into three com- 
ponents - the reproduction of the input state, the linear 
cross-talk and the circular-cross-talk. The measured out- 
put separated in this manner is shown in Figure [14l Only 
input U shows significant wavelength-dependent cross- 
talk. The overall degree of polarization, p, calculated as 
p 2 =q 2 + u 2 + v 2 , is very high for both Q and U inputs. 
The values range from 96-100% and thus the 'depolar- 
ization' terms are small. 

The negative and positive input states give results that 
are opposite in sign but very similar in magnitude and 
morphology. The linear cross-talk of the middle panel 
shows symmetry about zero and a very small 0.1-0.2 am- 
plitude. There is significant linear cross-talk for input U 
and a small amount for input Q at longer wavelengths. 
The third panel shows high UV cross-talk for input U 
at both short and long wavelengths with essentially a 
linear dependence with wavelength. The UV cross-talk 
fortuitously goes to zero around 6500A. 

The induced polarization is measured to be quite small, 
being under 5% to 6% in these wavelengths. Since the 
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induced polarization and depolarization terms are small 
with nearly 100% reproduction of the input degree of 
polarization, the Mueller matrix for the HiVIS system is 
essentially a 3x3 matrix involving QUV terms in rotation 
and linear-circular cross-talk. 

To investigate the image rotator further, the polariza- 
tion properties of the spectrograph are measured with 
the image rotator at angles of -90, -45, 0,-1-45 and +90 de- 
grees. Induced polarization (unpolarized flat field mea- 
surements) as well as purely polarized inputs were mea- 
sured at four polarizer orientations (+Q,+U,—Q,—U) 
for each of the 5 image rotator orientations. This set of 
20 observations allows us to characterize the wavelength- 
dependent changes with rotator angle. 

The first thing to note about the flat field is that there 
is significant linear polarization for all the image rota- 
tor orientations. The q terms all generally decrease to 
longer wavelengths while the u terms are more gener- 
ally flat with wavelength. The ±90° observations show 
a very significant u of 0.1-0.15 though with low v. The 
±45° observations show significant v at short and long 
wavelengths, with all observations showing zero induced 
circular polarization around 6200A, 300A shorter than 
other orientations. This shows that the induced polar- 
ization is at the 2% to 15% level and is highly dependent 
on the image rotator. 

The overall polarization properties of Figure [14] are 
generally reproduced at the varying image rotator ori- 
entations. The polarization reference frame rotates with 
the image rotator as expected. There is some change in 
the exact form of the linear to circular cross-talk but the 
relative amplitude does not change by much. The linear 
cross-talk is roughly double and there is some additional 
QV circular cross-talk when the rotator is oriented ±45°. 
We should note that some calibrations were repeated 
with the polarizer at different spatial locations and the 
resulting observations did not change significantly. 

5.2. Polarization Without Image Rotator 

Given the strong polarization variation when using the 
image rotator as well as the high fold angles of the mir- 
rors, we decided to implement a mode with the image 
rotator removed from the beam. This removes our con- 
trol of the projected slit orientation on the sky. However, 
we gain some very positive benefits. The removal of mov- 
ing oblique-fold mirrors greatly reduces the complexity, 
magnitude and wavelength-depenence of the cross-talk. 
We also gain in throughput and get reduced static wave- 
front errors caused by the high-angle folds. The AEOS 
telescope itself has five 45° folds before reaching our spec- 
trograph, but now there are no high-angle folds in the 
spectrograph optics. Though the rotator adds roughly 
15cm to the optical path after the collimator, this is a 
small change given the 550cm path length from the colli- 
mator to the tip tilt pupil image. The flat field illumina- 
tion pattern did not significantly change and the imaged 
spectral orders are essentially identical with and without 
the image rotator. The alignment of the system is good 
enough that the image rotator can simply be removed 
and replaced as needed without changing any other mir- 
ror orientations. 

The polarization measurements showed that the fore- 
optics induced polarization is typically 5% to 6% percent 
depending on the illumination pattern. The flat field 
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Fig. 15. — This Figure shows the measured polarization proper- 
ties with the image rotator removed. Pure input Q is shown with a 
solid line and was measured twice. Pure input U is the dashed line. 
The top panel shows the reproduced linear polarization input (Qi n 
to Qout, Ui n to Uout)- The middle panel shows linear cross-talk. 
The bottom panel shows the measured Stokes V running from to 
about 0.4. 

screen and tip-tilt guiding piezo were moved to mea- 
sure the response to a change in the illumination pat- 
tern. There are differences caused by the varying illumi- 
nation, but a coherent pattern is easily recognized. The 
flat field is dominated by Stokes q being around 0.05 to 
0.06. Stokes u is 5 to 10 times smaller, being 0.006 to 
0.012. The induced circular polarization, v, is another 
order of magnitude smaller than this. The total degree 
of polarization is essentially the degree of linear polariza- 
tion of 5% to 6% and this is essentially all Stokes q. This 
is quite a small variation given the 2% to 15% that can be 
induced with the image rotator at various orientations. 

The major improvements can be seen in the reproduc- 
tion of pure input linear polarization. Figure [15] shows 
the measured polarization properties for input Q (re- 
peated twice, solid lines) and U (dashed lines). The top 
panel shows the reproduced linear polarization input at 
a roughly constant 98% for input Q and a 98% falling 
to 90% for an input U. The middle panel shows the lin- 
ear cross-talk. The chromatic variation is very small and 
corresponds to less than 1°. The bottom panel shows the 
measured Stokes V which is nearly for input Q (just as 
with the image rotator in place) but only rising to 0.4 
for an input U (as opposed to ±1 changing to +1 with 
the image rotator in place). The general depolarization 
is much improved. This can be seen as the pure linear 
inputs being reproduced at the 90% to 98% level. The 
total degree of polarization is always above the 98% level. 
The mirror-induced rotation of the plane of polarization 
is nearly negligible, in contrast to the 10° seen with the 
image rotator. The induced V was very strong and very 
chromatic with the image rotator, but the induced V now 
has a maximum of 0.4 only for input U and only at long 
wavelengths. There is essentially no circular cross-talk 
at 5500A. 

As we had seen earlier, both the induced polarization 
and the depolarization are small. This alone suggests 
that one can assume that any linear to circular cross- 
talk has a corresponding circular to linear component. 
In order to verify this assumption, we remounted the slit 
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calibration optics (high-contrast polarizer and quarter- 
wave plate) just down stream of the diffusion screen in 
order to create a Stokes V input. We obtained several 
full-Stokes observations with both +V and — V inputs at 
various combinations of polarizer and quarter- wave plate 
orientations. The input V had U cross-talk in direct 
response to the UV cross-talk seen in Figure [151 We also 
measured linear polarization as a consistency check and 
found excellent agreement. 

As expected, the image rotator is the most signifi- 
cant source of induced polarization and cross-talk in our 
polarimeter. Removing the image rotator created sev- 
eral polarimetric advantages. The circular polarization 
cross-talk in the system is ±1 linearly changing to =pl 
from 5500A to 7000 A for linear polarization input at 45°. 
This induced circular polarization becomes much smaller 
without the image rotator, showing only UV cross-talk 
rising to 0.4 at long wavelengths. The chromatic rotation 
of the plane of polarization was up to 10° with the rota- 
tor but is less than 1° without the rotator. In general the 
depolarization is low both with and without the image 
rotator. The induced polarization without the rotator 
is almost entirely Stokes q at 5% to 6% with negligible 
induced circular polarization and with only a moderate 
variation with changing illumination patterns. This is 
to be expected as Stokes Q is defined parallel to the slit 
which is parallel to the image rotator axis. The spec- 
trograph is more easily deprojected without the rotator 
as the condition number of the deprojection matrix is 
much smaller. This allows us to recover an accurate po- 
larimetric reference frame at the entrance to the spectro- 
graph and to separate the telescope polarization proper- 
ties from the spectrograph. 

6. AEOS TELESCOPE POLARIZATION 

Knowing the polarization properties of the HiVIS spec- 
trograph with and without the image rotator, we can now 
quantify the polarization properties of the AEOS coude 
path. A simple experiment to illustrate the cross-talk 
is to look at twilight at the zenith. The linear polar- 
ization properties of twilight observed with HiVIS were 
investigated in Harrington & Kuhn 2008. They adopt a 
simple Rayleigh scattering model (cf. Liu & Voss 1997, 
Lee 1998, Cronin et al. 2006). With the Sun on the hori- 
zon in the west the zenith is highly polarized with the 
polarization angle pointing north-south. This is a simple 
approximation that holds quite generally (cf. Pomozi et 
al. 2001, Suhai & Horvath 2004, Horvath et al. 2002). 

The AEOS coude path has five 45° reflections and 
there is a pair of mirrors corresponding to the azimuth 
axis and a pair of mirrors corresponding to the altitude 
axis. If we point the telescope at the zenith and simply 
change azimuth through 360°, we can measure the polar- 
ization properties of this coude path as two mirrors cross 
and uncross. In Harrington & Kuhn 2008, the linear po- 
larization of the zenith was found to decrease from 70% 
to 80% at some pointings to 15% at other pointings. We 
now know that without the image rotator there is nearly 
no linear-circular cross-talk around 5500A but there is 
0.4 at 7000A. 

We repeated the twilight experiments without the im- 
age rotator using the new full-Stokes capabilities with 
achromatic wave plates. Figure [16] shows the measured 
Stokes parameters of twilight at the zenith on September 



5th 2009. On this night the Sun set at 280° azimuth. The 
left hand panel shows the Stokes parameters at azimuths 
of north, east, south and west (10, 100, 190 and 280). 
The right panel shows north-east, south-east, south-west 
and north-west (55, 145, 235 and 325). One can see that 
for the cardinal direction pointings of the left panel, the 
linear polarization is almost all Stokes q with v vary- 
ing from -0.2 to +0.2 from 5500A to 7000A. Conversely, 
at the non-cardinal pointings of the right panel, the ob- 
served polarization is substantially more chromatic and 
is entirely Stokes v at some wavelengths. The degree of 
polarization is again observed to be always 70% to 80% 
for all wavelengths. With these new observations we can 
conclude that the loss of linear polarization and the ro- 
tation of the plane of polarization seen in Harrington & 
Kuhn at 6500A is due to linear-circular cross talk. We 
are in the process of measuring the telescope polarization 
properties at all pointings. 

7. DISCUSSION AND CONCLUSIONS 

We have presented a Stokes-based deprojection 
method applied under various circumstances to highly 
chromatic astronomical systems. These methods are an 
easy and robust way to restore a reasonable polarization 
reference frame to a system with very significant chroma- 
tism and cross-talk. Since it is often more cost-effective 
to add spectropolarimetric capabilities to existing instru- 
ments, one typically has to deal with substantial system- 
atic effects. Stokes-based deprojection is a straightfor- 
ward and repeatable method to remove some of these 
errors. 

A simple laboratory spectropolarimeter was built to 
characterize and test our new optics as well as to demon- 
strate the Stokes-based deprojection methods in the lab. 
Retardances of our achromatic wave plates as well as 
LCVRs were measured. We verified the achromatic per- 
formance of the wave plates over their designed wave- 
length range. The LCVRs were characterized by creating 
a grid of retardances versus wavelength and voltage. An 
LCVR spectropolarimeter was constructed and a stan- 
dard sequence of observations were taken. The measure- 
ments were first deprojected using standard retardance- 
fitting. This retardance based deprojection was then 
compared to the Stokes-based deprojection method. The 
Stokes-based method was both easier to implement, less 
chromatic and more accurate. 

The Stokes-based deprojection methods were then 
demonstrated with a new HiVIS LCVR spectropolarime- 
ter. HiVIS now includes several optical / mechanical im- 
provements, a polarization calibration stage, a new full- 
Stokes achromatic wave plate mode and a liquid-crystal 
mode. This LCVR mode is now fully characterized to al- 
low the choice of polarization sensitivities optimized for 
a wide range of applications. We have shown that one 
need not map individual Stokes parameters to individual 
exposures. For instance, we are implementing a mode 
where our new fast shuffling detector can interleave or- 
ders on the focal plane. With this information, we can 
implement a mode where the condition number of the 
deprojection matrix is optimized for different sensitivi- 
ties within each exposure. Linear, circular or full-Stokes 
sensitivity can be specified within each frame to mini- 
mize systematic error in measurements of each Stokes 
parameter. 
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Fig. 16. — This Figure shows the measured twilight polarization properties using HiVIS without the image rotator. Each panel shows 
Stokes quv for twilight pointing at the zenith. The left panel shows cardinal pointings: north, east, south and west. The right panel shows 
non-cardinal pointings: north-east, south-east, south-west, north-west. The degree of polarization for all measurements was always 70% to 
80%. At non-cardinal pointings, the linear-circular cross-talk was severe and reached 100% at 6500A. 



We then used the achromatic wave-plate mode to fully 
calibrate the HiVIS spectrograph and show, as expected, 
that the main polarimetric effect of the spectrograph mir- 
rors is linear-circular cross-talk caused by the image rota- 
tor. With this image rotator removed, the spectrograph 
still has linear-circular cross talk at longer wavelengths at 
the 0.4 level. We then used this spectrograph calibration 
information to show that the AEOS telescope also in- 
duces 100% linear-circular cross-talk at the zenith when 
at non-cardinal azimuths. The image rotator and tele- 
scope can both induce 100% linear-circular cross-talk. 
These polarimetric results are repeatable and as such, 
may be calibrated. Accurate absolute polarimetry from 
any modern alt-az telescope requires careful cross-talk 
calibration as we describe here. 

Stokes-based deprojection applied to highly chromatic 



systems will significantly improve the repeatability and 
accuracy of spectropolarimetric measurements. The 
LCVR system, though nominally chromatic, can be used 
to obtain broad-band spectropolarimetry with good ac- 
curacy and repeatability. This achromatization allows 
one to realize the performance advantages of a fast- 
switching system with no moving parts. We outlined how 
systematic errors resulting from rotating wave-plates and 
slow modulation may now be eliminated from the HiVIS 
spectropolarimeter. 
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